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Trapped singly-charged ions can crystallize as a result of laser cooling. The emerging structure 
depends on the number of particles and on the geometry of the trapping potential. In linear multipole 
radiofrequency traps, the geometry of the radial potential can lead to the formation of single-ring 
structures. We analyse the conditions and stability of single rings as a function of the number of 
poles. The rings form tubes for a large number of ions and sufficiently small trap aspect ratios. 
In these structures the arrangement of the ions corresponds to a triangular lattice folded onto a 
cylinder. The stability of the tubular structures is numerically studied for different lattice constants 
and their normal mode spectrum is determined. 



I. INTRODUCTION 

Singly-charged particles in an external trapping poten- 
tial can form ordered structures at sufficiently low tem- 
peratures, for which the interaction energy exceeds the 
thermal energy [U HJ. The structures emerge from the 
balance between the Coulomb repulsion and the confining 
force, and their form can be thus controlled by changing 
the geometry of the trap [JJ. Long-range order in struc- 
tures of thousands to millions of ions were experimentally 
demonstrated in Penning and linear Paul traps [3J H] . 

While most works have focussed on harmonic traps, 
or better said, potentials which can be approximated 
by harmonic ones, most recently research has focussed 
on dynamics and crystallization in anharmonic poten- 
tials formed by multipole traps |S]. Recent works have 
studied the formation of ordered structures [5llll| and 
crystallization in octupole traps has been observed by 
means of laser cooling [7] . These studies show that ions 
in multipole traps tend to form hollow structures at very 
low temperatures. In linear multipole traps the struc- 
tures take the form of either rings, tubes, or cylinders 
formed by several layers of tubes. The degree of ordering 
in these configurations depends on the trap parameters 
and on the number of ions: Numerical simulations indi- 
cate that both chirality and localized defects can appear 
as a result of a mismatch between the number of ions and 
the "magic number" that would correspond to a perfect 
pattern for the given trap parameters [TU]. Applications 
are being discussed for metrology and quantum simula- 
tors [9]. This perspective further motivates a detailed 
characterization of these systems. 

In this work we theoretically analyze the conditions 
under which single rings form in a multipole linear trap, 
like the one sketched in Fig. [TJ and study their stability 
as a function of the number of poles. The conditions for 
stability of a single, two, or multiple rings are determined. 
The study is extended to a large number of ions for the 
purpose of identifying the parameter regime for which 
ion tubes, namely, an ordered array of single rings, are 
stable. The eigenmodes and eigenfrequencies of the tubes 
are then evaluated. 

The article is organized as follows: In Section ITT] we in- 



troduce the effective potential for linear multipole traps. 



Section III is devoted to the analysis of the stability of 
single rings in such potentials. In Section |IV| the con- 
ditions for the stability of tubular structures are inves- 
tigated and the normal modes are determined. Finally, 
in Section |V] we summarize the main results and discuss 
some outlooks to this work. The appendix provides fur- 
ther details relevant to the discussion in Sec. Mil 



II. THE MODEL 

In this section we introduce the basic Hamiltonian that 
determines the dynamics of N singly-charged particles 
confined by a linear multipole trap. Some preliminary 
considerations on the charge distribution in presence of 
an anharmonic potential are made. 



A. Hamiltonian 

The particles are all assumed to possess the same 
charge q and mass m. They mutually repel via the 
Coulomb interaction and are subject to a trapping po- 
tential which we denote by Vt rap (r) and is a function of 
the position r. At sufficiently low density the particles 
are distinguishable and we label them by j = 1, . . . , JV, 
so that Tj and Pj are the canonically-conjugated posi- 
tion and momentum of ion j. The energy of the charges 
is given by Hamiltonian 



^ Pi 



.i=i 



2m 



N 



N N 



1 



47re f-f , 4? 



= 1 l=j+l 



ri 



with eo the vacuum permittivity. 

The potential Vt rap (r) is generated by a radio- 
frequency trap. The trap consists of a set of 2k parallel 
electrodes that are equally spaced on a cylinder with ra- 
dius r , as shown in Fig. [I] for 2k — 8. To provide the 
confinement along the z-axis, the electrodes can be split 
into three segments, and a positive static voltage is ap- 
plied to the endcaps of the electrodes [6] . For small volt- 
ages and in the central region of the trap, the resulting 
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FIG. 1. Schematic representation of a linear octupole ra- 
diofrequency trap. The trap consists of eight parallel rods 
that are situated on a ring with equal radial spacing. Each 
rod is separated into three segments. An oscillating voltage 
is applied to the adjacent parts of the central segments of the 
electrodes and a static voltage to the outer segments of the 
electrodes [6]. 

static potential can be approximated by the equation 
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with uj z the axial trapping frequency, which depends on 
the static voltage applied, while r = y/x 2 +y 2 is the 
radial coordinate and z the position along the trap axis. 
The radial confinement is generated by a rapidly oscillat- 
ing voltage, which is applied to the central part of each 
electrode. Such setup creates a potential that in cylin- 
drical coordinates takes the form [5J [T2] 



$ rf (r, 6) = V 1 — 1 cos(fc6») cos (fit) , (3) 

where f2 is the angular frequency of the oscillation, 6 is 
the azimuthal angle, t is the time, and Vq is a real pos- 
itive scalar with the dimensions of an energy, such that 
the amplitude difference between two adjacent electrodes 
is 2Vb- For k > 2 a general solution of the equations of 
motion is not known. Nevertheless, in some cases the 
orbits of a single ion in potential ^ can be well ap- 
proximated by the superposition of a slow motion along 
a smooth curve and a fast oscillation about it at much 
smaller amplitude, as it can be done in the stability re- 
gion of a Paul trap. This is verified provided that the 
relative variation of the local electric field, seen over the 
amplitude of the oscillations of the driven motion, is suf- 
ficiently small, namely, when rj < 0.3, where rj is the 
so-called adiabaticity parameter defined as 
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and E is the local field generated by the oscillating po- 
tential [13l [14] . In this limit the secular motion is subject 
to an effective conservative force associated to the elec- 
trostatic potential usually denoted by pseudopotential, 



which in this case reads [T3] 
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with v}f = q Vo fc/(mO Tq). The effective trap potential, 
which confines a particle in all three directions of space, 
is then given by 



V tIap (r)=V li (r) + V st (r), 



(6) 



and is the potential we consider in Eq. ([I]). It neglects 
any influence of mirror charges on the conducting bound- 
aries of the electrodes, which is a good description if the 
ions are sufficiently far away from the electrodes pQ. 

Potential Vtrap(?") can be written as a sum of two 
contributions, one depending on the radial coordinate 
r and the other on the axial coordinate z. While in a 
quadrupole trap (fc=2) the pseudopotential is harmonic 
and its minimum is at the trap center (which here coin- 
cides with the origin of the coordinates), in a multiple 
trap, for k > 2, the center of the trap is unstable and the 
potential is minimum along a cylinder centered at the 
trap axis with radius [51 Q3] 
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When more than one ion is inside the trap, the Coulomb 
interaction must be taken into account to determine the 
equilibrium density. 



Charge distribution at equilibrium for weak 
correlations 



The distribution of charges at equilibrium in a mul- 
tipole trap has been numerically evaluated in a series 
of works, which have reported the appearance of hollow 
structures [SHU) . This phenomenon can already be un- 
derstood by calculating the charge density distribution 
of a weakly correlated plasma |14j . For particle numbers 
of the order of hundreds and more, the thermal equi- 
librium states can be described by a Gibbs distribution. 
The particle density, which we denote by n(r), can be 
approximated by the single-particle Boltzmann distribu- 
tion in the limit in which the Debye length A# is small 
with respect to the plasma dimensions [T|. The density 
takes the form 



n(r) ~ Aexp{-f3[V tlap (r) +q^ p (r)}}. 
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where (3 = 1/ksT, with ks Boltzmann constant and T 
the temperature, A is a scaling constant warranting that 
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N, and is the electric potential created 



by the charge distribution itself, satisfying the Poisson 
equation 



A$ p (r) = -—n(r) . 
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Potential $ p is assumed to vanish at the conducting sur- 
faces of the trap, and infinitely far away from the charge 
distribution (these boundary conditions allow for neglect- 
ing the mirror charges only when the conducting surfaces 
are taken to be far away from the plasma) . When the po- 
tential is harmonic the density is uniform, n(r) = hq, and 
the thermodynamic properties are solely characterized by 
the so-called plasma parameter T — q 2 n^ 3 / (AneokBT). 
The Debye length takes the form A,d = ^huTe®/ {q 2 n$), 
and is larger than the interparticle distance when Y < 1, 
since n X 3 D = (47iT;r 3 / 2 . 

In multipole traps the scenario is quite different. Solv- 
ing Eq. {9) in presence of an anharmonic potential leads 
to non-uniform densities: the density hrst increases with 
the distance to the trap axis, reaches a maximum at a 
non-zero radius, and then rapidly drops to zero. We refer 
the reader to the appendix and Ref. [H] for detailed cal- 
culations. This behaviour reflects the form of the radial 
dependence of the potential for multipole traps. One can 
dehne a plasma parameter T(r) = o 2 n(r) 1 / 3 / (47reo^s7 1 ), 
which results to be position dependent, and is hence 
smaller at the center of the trap, reaching a maximum 
at the radius where the density is maximum. This fea- 
ture indicates that crystallization is first expected in the 
outer shells of the crystals, and only at sufficiently low 
temperatures the ions form ordered clusters in the entire 
trap volume where they are distributed [HI [TO] • 



Charge distribution at equilibrium for strong 
correlations 



For uniform distributions in three dimensions a singly- 
charged plasma exhibits at V ~ 174 a first-order phase 
transition from a liquid to a crystal phase [h . Similar to 
the case of harmonic traps, in anharmonic traps at suf- 
ficiently low densities and temperatures the ions crystal- 
lize. By this, we mean they arrange in structures in which 
each ion has a well-defined equilibrium position and 
performs small oscillations about it, whose amplitudes 
are much smaller than the interparticle distance. In an 
equilibrium configuration the net force over each ion van- 
ishes, namely, 



dV 
Or* 



= 0, 



where V is the total potential, corresponding to the po- 
tential term in Eq. 0, which we here rewrite as 
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The set of positions {r^} describes a stable configura- 
tion if the potential there evaluated is a local minimum, 
which corresponds to the condition that the Hessian of 
the potential energy at the positions rj ' is positive- 
definite. We note that, in general, there will be several 
possible local minima for a given number of system pa- 
rameters. Unless otherwise stated, we will discuss the 
configuration that corresponds to a global minimum of 
the potential energy. 



III. ION RINGS 

In this section we analyze the stability of single rings 
which, for a given number of ions, are formed at the cen- 
ter of the trap for a sufficiently large aspect ratio between 
the radial and the axial trapping potential. 

In order to study the problem, we first introduce di- 
mensionless coordinates, and denote by r = r/f the di- 
mensionless position, with f = [axj 'a3)~ 1 ^ 2k ~ 1 ' . The 
dimensionless potential V is obtained by rescaling po- 
tential V in Eq. (10) by the quantity 0,3 /f and reads 
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where the parameter 

c = a 2 f 3 /a 3 (13) 

depends on the relative strength of each contribution. 
Equation (12 1 shows that every cluster is determined by 



three independent parameters, which can be chosen to 
be N, k, and c. In the rest of this section we discuss the 
equilibrium configurations and omit the tilde (but always 
refer to the dimensionless coordinates, unless otherwise 
stated). 

The equilibrium configurations are found by numeri- 
cally minimizing the potential energy in Eq. ( 12 1 using 



(11) 



simulated annealing 15, Hi . Before we start, we remark 
that single ring structures in linear multipole traps have 
been reported in Rcfs. [6l[9l[TT]. In Ref. [9] the condi- 
tions for stability of a single ring were derived and in Ref. 
[11] numerical simulations based on molecular dynamics 
were reported. The results shown below are in agreement 
with these studies and complement them since they pro- 
vide a systematic characterization of the transition from 
a single to a double ring by means of a closed set of equa- 
tions. 



A. Single rings 

We first focus on the existence of single rings in the 
plane z = of the multipole trap. Given a number of 
ions N in a multipole trap with k > 2, one can always 
find a sufficiently large value of c for which the ions ar- 
range on a single ring at the plane 2 = 0, where the ring 
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FIG. 2. Equilibrium configurations of iV = 12 ions in an 
octupole (k = 4) trap, determined by simulated annealing. 
The panels on the left show the charge distribution in the 
radial, x — y, plane for different values of c, the panels on the 
right the corresponding distribution on the plane x — z. In 
(a)-(b) the value c = 10 was taken. The plots in (c)-(d) are 
evaluated for c = 8.5 and exhibit a double ring structure. 



radius R depends on the three parameters N, k, c. Some 
configurations are shown in Fig. [2] Analytical expres- 
sions for the radius have been derived in [9] and will be 
reported later on for studying the transition to the double 
ring visible in Figs. j2jc) and (d). 

Few remarks are here in order. For sufficiently large 
aspect ratios between the radial and the axial potentials 
(corresponding to large values of c) , in a linear Paul trap 
single-ring structures are a global minimum of the poten- 
tial only for 1 < N < 5 [17] ■ In multipole linear traps, 
instead, these configurations are found for any number 
N provided that c is sufficiently large. They result from 
the specific shape of the anharmonic potential in Eq. ^ , 
having a minimum shifted away from the trap center. 
Indeed, even for a single ion the radius r m j n , Eq. Q, 
increases as c is increased either by decreasing the ra- 
dial confinement (and thus a\) or by rendering the axial 
confinement steeper (increasing 02). 



B. From a single to a double ring 

We now analyze the ions' configuration as c is de- 
creased keeping k and N constant. In this case, a single 
ring in the plane z = is stable until c reaches a crit- 
ical value c cr it, below which the ions' distribution is no 



FIG. 3. Equilibrium configuration of 17 ions in a linear mul- 
tipole trap for k = 4 and c = 10, determined by simulated 
annealing. The ions form two rings with equal number M = 8 
of ions, while one ion is located in the plane z — between 
both rings, forming a localized defect. 
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FIG. 4. Equilibrium configuration of 70 ions in an octupole 
trap (k — 4) for c = 4, corresponding to a tubular arrange- 
ment. The positions have been found by simulated annealing. 
The inner rings are composed by 12 ions per ring, while at 
the two edges the rings are composed by 10 ions each. 



longer planar. An analogous behaviour is observed if c 
is kept fixed but the number of ions N is increased: in 
this case the transition to a three-dimensional structure 
is observed when exceeds a value A cr i t which depends 
on k and c. 

Be now N constant and c below (but sufficiently close 
to) c cr it: For an even total number of ions, N = 2M, 
the ions arrange in two rings with the same number of 
ions and the same radius, as shown in Fig. [2]jc)-(d). In 
Fig. |2jd) one observes that the structure seen from the 
side takes a zigzag shape. This kind of shape was also 
reported in Refs. HI]- For an odd number of ions, 
N = 2M + 1, the ions form two rings with equal number 
M of ions, while one ion sits in the plane z — between 
both rings, forming a localized defect, as shown in Fig. [3] 

By further decreasing c, parameter regimes are found 
for which the ions arrange in multiple rings forming a 
cylinder, as displayed in Fig. [4] In particular, all rings 
have the same number of ions and approximately the 
same radius, except for the rings at the edges that have 
a smaller radius and also a smaller number of ions. 

In the following we derive an analytical expression for 
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the instability of a single ring composed by an even num- 
ber of ions. Close to the instability point, the following 
ansatz describes the ions' equilibrium positions either in 
the single or double ring configuration: 



Xj = i?cos(27T j/N) 
y,= Esin(27r j/N) 
Zj = (-If d/2 



3 = h 



,N, 



(14) 



where R is the (dimcnsionless) ring radius and d/2 is the 
(dimcnsionless) distance of the ions from the plane z = 0. 
Therefore, if d = 0, all ions are located in the plane z = 
forming a single ring of radius R, while for d ^ 0, the ions 
form two rings, each consisting of A/2 ions, with radius R 
and at distance d from each other. The potential energy 
of a double-ring configuration depends on the radius R 
and distance d between rings, and is obtained by using 
Eqs. (fill) in Eq. (pi 
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We now derive the relations that allow us to find the 
critical value c cr it for which the transition from a sin- 
gle ring to a double ring occurs. If a single ring with 
radius R = R$ is stable, then the potential energy will 
have a minimum at d = 0. This implies that the first 
derivatives of the potential energy with respect to R and 
d have to vanish, 
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where the second condition always holds for symmetry 
reasons. The second derivatives must satisfy: 
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while the mixed derivative d 2 V/(dRdd) vanishes at 
d = 0. 

At the point where a single ring becomes unstable and 
splits into a double ring, the second derivative with re- 
spect to d changes its sign crossing a zero, 



FIG. 5. Critical parameters i? C rit (a) and c cr i t (b) at the 
transition point from a single to a double-ring structure as a 
function of the particle number for an octupole trap (k = 4). 
The curves are evaluated from the analytical expressions in 
Eqs. JT8l)-(pT|. 
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>0, (16) 

The corresponding value of the ring radius reads 



^?crit 



2S 1 (N) + S 3 (N) 



16(fc- 1) 



l/(2fe-l) 



d 2 V 



dd 2 



0. 



(17) 



(Ro,0) 



To find the point where a single ring becomes unstable, 
we solve Eqs. (15) and (17) and find that the parame- 
ter Ccrit reads 



Figure [5] displays the parameters c cr ;t and i? cr it as a func- 
tion of the particle number N in an octupole trap (k = 4). 
For large particle numbers, the parameter c cr it scales as 
Ccrit ~ jv 6 Cfc+i)/(2fc-i) > while RcM „ jv3/(2fc-i)_ Thege 

behaviours show that in a quadrupole trap the critical 
radius scales linearly with the particle number, while for 
high orders of the multipole potential, k 1, the critical 
radius reaches a constant value. In fact, for large values 
of k, the multipole potential approaches a box potential 
of the form 



Ccrit = ^S 3 (N) 
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which has thus maximal radius Rq = 1 (in dimension- 
less coordinates). Note that Eq. (21) exhibits a differ- 



ent dependence on A than the critical radius derived in 
Rcf. [9 . This is due to the different conditions which 
were considered. In [9] the axial trap frequency is fixed. 
Here, we assume that the radial confinement is kept con- 
stant, while c is varied by changing the axial trap fre- 
quency: Correspondingly, Eq. (21) (once rescaled by 
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r = (a 3 /ai) 1 ^ 2k ^) is reported in terms of a\, 03, fc, 
and N. 



C. Discussion 

The transition from a single to a double ring shares sev- 
eral analogies with the linear-zigzag transition observed 
in a linear Paul trap [T8l420| . Indeed, as for the linear- 
zigzag transition a symmetry is broken, which is the sym- 
metry by reflection about the z = plane. However, in 
the case of a multipole trap the radius of the distribution 
changes as the aspect ratio varies across the structural 
transition, and thus the density along the direction de- 
fined by the polar angle varies. An expansion close to the 
critical point shows that the distance d increases like y/Sc 
for 8c — c cr u — c > 0, while R decreases linearly in 5c. 
In addition, the effect of the curvature of the ring is rel- 
evant as long as the radius remains finite. Nevertheless, 
this system can be considered a platform where to test 
predictions related to quantum effects at the structural 
transition [2TJ[22]. For this purpose, an adequate analysis 
following the lines of the study in [531 HI] is required. 

Figure [4] shows an example of configuration which is 
found by further lowering the value of c. For a very 
large number of ions and weak confinement along the z 
axis the ions form tubes, namely, uniformly spaced rings 
with the same radius and number of ions per ring. Their 
characterization is the subject of the next section. 

IV. ION TUBES 

In this section we study the stability and dynamics 
of structures composed of multiple rings forming tubular 
structures. These structures have been first predicted 
in Ref. 16] . The arrangement of the ions in the central 
region of such tubes corresponds to a triangular lattice 
folded onto a cylinder. We note that previous studies 
based on molecular dynamics simulations have observed 
the formation of inner rings for systems with numbers 
of ions of the order of one hundred [B] to one thousand 
|llj . Besides, the impossibility to fit any number of ions 
into a regular pattern for given trap parameters has been 
shown to lead to chiral structures or localized defects [10] . 
In the following we will restrict to the study of tubular 
structures with no inner rings. Furthermore, we focus on 
the configurations with neither chirality nor defects. This 
means that we consider large ion numbers with very weak 
axial trapping for which single tubes are formed [TO] . 

A. Equilibrium structure 

Let the tube be composed by 2Na rings, each with 
radius R and Nr ions per ring, and with distance d/2 
between adjacent rings. The tube is sketched in the up- 
per panel of Fig. [6] Periodic boundary conditions in z 



direction are assumed, taking that the trap confines the 
ions only in the radial direction, 

ot 

while the axial confinement vanishes, 0,2 = 0, so that the 
total potential V the ions feel reads 

trap 4^eo 4l r a -r a , 

This assumption can be taken for a long tube in a linear 
multipole trap for which the boundary effects can be ne- 
glected. Under this assumption, the resulting structure 
possesses a discrete rotational symmetry by the angle 
4>o = 2tt/Nr about the z axis, as visible in the lower 
panel of Fig. [6j and a discrete translational symmetry 
under the axial displacement d. These symmetries lead 
to the choice of a unit cell composed by two ions, as 
shown in the lower panel of Fig. [6] The total number of 
unit cells is iV ce u s = Na Nr, thus half of the total number 
of ions composing the tube. 



1-1 I 1 + 1 




FIG. 6. Top panel: sketch of the ring array at uniform 
distance d/2 along the z axis. The lower panel shows the ions' 
distribution in two adjacent rings: this structure is repeated 
along the chain. An example of the chosen elementary cell 
is given by the pair of ions labeled with a. The other labels 
are j = 0, . . . , Nr — 1 for the number of cells in two adjacent 
rings and 1 = 0,..., Na — 1 for the pair of rings. 

In what follows we shall label the ions by the set of 
indices a — (j, I, a): the index j labels the cells along the 
azimuthal angle with j = 0, . . . , Nr — 1, the index I labels 
the cells along the z direction with I = 0, . . . , Na — 1, 
while a — 0, 1 labels the two ions in each unit cell, see 
also Fig. [6] The equilibrium positions of the ions are 
given by the vectors R a = R(j,l,cr), which in cartesian 
coordinates read 

/i?cos [fo (j+a/2)}\ 
R a =\Rsm[<j>o(j + a/2)] \. (24) 
\ d{l + a/2) J 
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Here, R is found from the condition that the net force on 
each ion vanishes, and fulfills the equation 
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Note that the sum on the right-hand side converges for 
an infinitely long tube. 

If this configuration is stable, at sufficiently low tem- 
peratures each ion performs small oscillations with a dis- 
placement u a from its equilibrium position R a , so that 



Ra 



(26) 



where r a indicates the position of the particle labelled by 
a and where periodic boundary conditions impose that 
u jd,a = Uj+N R ,l+N A ,o- (with u a = Uj t i <(T ). The harmonic 
approximation can be performed when the displacement 
is at all times much smaller than the interparticle dis- 
tance. In this limit the potential energy can be approxi- 
mated by its quadratic expansion 



2 ^ 



u a D(a, a')u a 
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TABLE I. Number of separate stability regions of the tubes 
for different multipolar orders k and number of ions per 
ring Nr, taking Nr > 4. For low Nr two stability regions 
are found. 
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For Nr > 14 no stable region was found. 



finds one or no stability region. The number of different 
stability regions is given in Table [i] as a function of the 
multipole order k. 
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where Eq — V'({R a }) is the energy of the equilibrium 
configuration, while D is a matrix whose coefficients read 



D^ v (a,a') 



d 2 V 



G {x,y,z} 
(28) 



and all derivatives are evaluated at the equilibrium posi- 
tions {R a }- 

The analysis of the eigenvalues of matrix D as a func- 
tion of the axial distance d and of the trap parameter a\ 
allows us to determine a stability diagram for the tubu- 
lar structures. We find the stability regions by identi- 
fying the parameters for which the eigenvalues of D are 
positive. The borders correspond to the case where one 
(or more) eigenvalue vanishes, while outside of the sta- 
bility region there are negative eigenvalues, i.e., the tube 
is unstable. The stability diagrams were calculated for 
up to 20 ions per ring for different multipolar orders k. 
Figure [7] displays one such diagram for the case of an 00- 
tupole trap. Depending on d and a\ there are different 
numbers of distinct regions where the tubes are stable. 
The corresponding stability boundaries are linear when 
plotted as a function of the radius R and of the distance 
d. This behaviour can be understood in terms of scal- 
ing arguments: in fact, the equilibrium configuration can 
be rescaled by a factor A if the trapping parameter a\ is 
rescaled by a factor A~ < - 2fe ~ 1 '. In general, for smaller Nr 
and higher multipole orders k there are up to two differ- 
ent stability regions, while for large Nr and small k one 



FIG. 7. Stability diagram of ion tubes with six ions per ring 
(Nr — 6) in an octupole trap (k = 4), here in the parameter 
space d and a\. The shaded areas indicate the regions where 
the ion tubes are stable. The value of tube radius depends on 
both d and ai through Eq. | |25[ ). The inset displays a zoom 
where the parameters of the plots in Fig. [8] are shown. 



B. Normal modes 

In the rest of this section we discuss the normal modes 
of the ion tube. For this purpose, we employ a coordinate 
system that matches the symmetry of the crystal. Thus, 
we choose for each particle a a coordinate system such 
that x a denotes the coordinate in the radial direction 
and y a the coordinate in the azimuthal direction. The 
potential energy in the new coordinates reads |25j 



V = E + i ^ u a D(a, a')u a i 

ct.a' 



(29) 



where 



Ua = Q~ 1 (a)u a (30) 
are the coordinates in the new rotated system and 

b(a,a')=Q- 1 (a)D(a,a')Q(a') (31) 



are the new matrix coefficients. The rotation matrix used 
here is given by 



Q0» = 



fcos(4>ja) -sin(^>) 0" 
sin(^ CT ) cos^^) 





(32) 



where 



f>o(j + c/2). Clearly, the transformation 



does not depend on the index I labeling the pair cell 
along the z axis. The equations of motion for each ion 
read 



mu 



i a = — D(a, a!)u a i. 



(33) 



Because of the periodic boundary conditions, this equa- 
tion can be solved using the ansatz 



a(*i.*».») = eo .exp 



2wkij 27rfc 2 i 

1 I 1 ^ w fci,fe 2 ,n^ 



N A 



(34) 



where k\ and &2 take the discrete values k\ — 0, . . . , Nr — 
1 and k 2 = 0, . . . , N A -1. Plugging Eq. p4h into Eq. Q 



leads to an eigenvalue problem, such that for every pair 
of values (fci,^) one finds the frequencies u>ki,k2,n> with 
n = 1, . . . , 6, whose corresponding eigenvectors are e^. 

The eigenfrequencies were calculated numerically for 
different values of A^, Nr, d, and R. Typical spectra 
for a fixed distance d = 10/xm are shown in Fig. [8j The 
spectra are evaluated for 400 rings, but we note that, by 
comparing the spectra for different values of the number 
of rings Na, one finds that tubes composed by tens of 
rings have a spectrum which is already close to the one 
reported here. 

Each row in Fig. [8] corresponds to a different steep- 
ness of the radial trapping potential, a±, while each col- 
umn corresponds to the spectra obtained for values of 
k\ = 0,±1,±2,±3. In particular a\ decreases from top 
to bottom: the value in the first row (from above) corre- 
sponds to the parameters of the upper point in the inset 
of Fig. [7J the second row corresponds to the parameters 
of the central point, and the third row to the parameters 
of the point at the lower edge of the stability region. 

We observe the existence of 6 branches for each value 
of k\. They are reduced to three (each twofold degener- 
ate) for k\ = ±3. Because of the continuous rotational 
and translational symmetries of the trapping potential, 
two eigenfrequencies at fej = k-z = are zero: one of 
these corresponds to the rigid rotation about the z-axis 
of the crystal, the other to rigid translations along z. In 
general, the eigenmodes at k\ = can be understood in 
terms of normal modes of a two-dimensional biperiodic 
crystalline structure For fc 2 = 0, on the other hand, 
the eigenmodes are independent of the axial coordinate 
z. We denote the eigenmodes with k\ = fc 2 = by bulk 
excitations. 

The tube becomes unstable at the points where eigen- 
frequencies, other than the bulk eigenfrequencies, vanish. 



For the cases we studied the corresponding modes do not 
offer an intuitive interpretation. For some multipole or- 
ders and numbers of ions per ring, they correspond to 
motion in azimuthal and axial direction suggesting a re- 
ordering of the structure within a cylinder without chang- 
ing the radius. 



V. CONCLUSIONS 

The stability of ion rings in linear multipole traps was 
studied both analytically and numerically. Two types of 
clusters were considered: single rings, which can form in 
the center of the linear multipole trap, and ion tubes, 
which arc found for large numbers of ions when the ax- 
ial potential is sufficiently flat. For the latter case, the 
stability diagram and the normal mode spectrum in the 
stability region were determined. 

This study sets the basis for a thermodynamic char- 
acterization of these structures in the spirit of the work 
in Ref . [2 3) . Such study will require the definition of an 
appropriate thermodynamic limit, which can be identi- 
fied with the number of rings going to infinity, Na — > oo, 
while the distance d between rings and the radius R re- 
main constant, meaning that the linear density remains 
constant. A different thermodynamic limit can be identi- 
fied for a single ring structure. This will in fact approach 
a linear chain in the theoretical limit in which the radius 
of the ring becomes infinite while the ions density on the 
ring is kept constant. An experimental situation deals 
with finite size systems, which include the effect of the 
curvature of the ring. For odd numbers of ions, the insta- 
bility of the ring leads to the creation of localized defects, 
which can be considered for quantum information plat- 
form with localized kinks, in the spirit of the proposal in 
Ref. [H]. 

The influence of quantum fluctuations at the instabil- 
ity between a single and a double ring, furthermore, can 
provide additional aspects to the one identified in Ref. 
[2 2) for the case of a linear chain in a linear quadrupole 
trap. 
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FIG. 8. Spectra of 400 rings of Ca + ions in an octupole trap for different trap parameters. Each ring is composed by six 
ions, the distance between rings is 10/im. The first row, (a)-(d), shows the spectrum for oi = 3.4 • 10 9 J m -6 (upper point in 
the inset of Fig. [7]), the second row, (e)-(h) for a 1 — 8.4 • 10 8 J m~ 6 (central point in Fig. [7]), and the third row (i)-(l) for 
ai = 2.7 • 10 s J m -6 (lower point in Fig. m). Each column corresponds to a value of fci: from left to right fci = 0, ±1, ±2, ±3. 
Because of the specific symmetry of the crystal, for fci = Nr/2 (when Nr is even) each frequency is twofold degenerate at a 
fixed value of ki. 



Appendix A 

In order to find the density profile, the two coupled 
equations ([8|-(|9| need to be solved. This can be done 
in a straightforward way in the case of an infinitely long 
trap |14j . It is convenient to introduce the logarithmic 
density profile 



*(r) = In 



n(r) 
ra(0) 



(Al) 



and express the density as a function of this quantity, 

n(r) = 71(0) exp [*(r)] . (A2) 
By inserting eq. ( A2 ) in the Poisson eq. ([9| , one finds 
1 I 



A*(r) 



A 



(0) 



exp [*(r)] 



q 2 n(0) 



AV tIS 



(A3) 



where we call \p the Debye length in the center of the 
trap. We note that this equation only depends on the ef- 



fective potential of the trap through the Laplacian; since 
the static potential V s t must satisfy the Laplace equation, 
it has no influence on the density profile within the region 
of validity of approximation (fSl) . In rescaled cylindrical 



coordinates p 
explicit form of Vrf 



r/\$ and g 



z/Ajp , and writing the 



Eq. (A3 1 becomes 



1 d d 
p dp^ dp 



#(p,0=exp[¥(p,0]-7P 



2fe-4 



with the parameter 

= ^ — 
no 4m 



Vk{2k-2) 



fir* 



v (0) 



2fe-4 



(A4) 



(A5) 



For a very long cylinder, the £ dependence can be ne- 
glected and the equation can be integrated. The results 
of this integration are shown in Fig. [9j While for small 
Debye lengths (small 7) in a quadrupole trap (2k = 4) 



10 




(a) 7 = 0.029, 2k = 10 (b) 7 = 1.5, 2k = 4 



FIG. 9. Density profiles in (a) a multipole trap with k — 5 and (b) a quadrupole trap (k = 2) for different values of the 
parameter 7. 



there is a nearly constant density in the center of the trap, 
in any multipole trap (2k > 4) the density increases from 
the center, reaches a maximum and then drops to nearly 
zero within one Debye length. Large Debye lengths, cor- 



responding to relatively high temperatures of the plasma, 
result in density profiles in multipole traps which are 
nearly constant over the whole charge distribution. 



[1] D. H. E. Dubin and T. M. O'Neil, Rev. Mod. Phys. 71, 
87 (1999). 

[2] E. Wigner, Phys. Rev. 46, 1002 (1934). 

[3] J. N. Tan, J. J. Bollinger, B. Jelenkovic, and D. J. 
Wineland, Phys. Rev. Lett. 75, 4198 (1995); T. B. 
Mitchell, J. J. Bollinger, D.H. E. Dubin, X.-P. Huang, 
W. M. Itano, and R. H. Baughman, Science 282, 1290 
(1998); W. M. Itano, J. J. Bollinger, J.N. Tan, B. Je- 
lenkovic, X.-P. Huang, and D. J. Wineland, Science 279, 
686 (1998). 

[4] M. Drewsen, C. Brodersen, L. Hornekaer, J. S. Hangst, 

and J. P. Schiffer, Phys. Rev. Lett. 81, 2878 (1998); 

A. Mortensen, E. Nielsen, T. Matthey and M. Drewsen, 

Phys. Rev. Lett. 96, 103001 (2006). 
[5] D. Gerlich, Adv. Chem. Phys. 82, 1 (1992). 
[6] K. Okada, K. Yasuda, T. Takayanagi, M. Wada, H. A. 

Schuessler, and S. Ohtani, Phys. Rev. A 75, 033409 

(2007). 

[7] K. Okada, T. Takayanagi, M. Wada, S. Ohtani, and H. 

A. Schuessler, Phys. Rev. A 80, 043405 (2009). 
[8] F. Calvo, C. Champenois, and E. Yurtsever, Phys. Rev. 

A 80, 063401 (2009). 
[9] C. Champenois, M. Marciante, J. Pedregosa-Gutierrez, 

M. Houssin, M. Knoop, and M. Kajita, Phys. Rev. A 81, 

043410 (2010). 

[10] E. Yurtsever, E. D. Onal, and F. Calvo, Phys. Rev. A 

83, 053427 (2011). 
[11] M. Marciante, C. Champenois, A. Calisti, and M. Knoop, 

Appl. Phys. B 107, 1117 (2012). 
[12] P. K. Ghosh, Ion traps, Clarendon Press (Oxford, 1995). 
[13] E. Teloy and D. Gerlich, Chem. Phys. 4, 417 (1974). 



[14] C Champenois, J. Phys. B: At. Mol. Opt. Phys. 42, 
154002 (2009). 

[15] L. Ingber, Mathematical and Computer Modelling 18, 29 
(1993). 

[16] T. F. Coleman and Y. Li, SIAM Journal on Optimization 
6, 418 (1996). 

[17] In quadrupole traps, for N > 5 the single-ring configu- 
ration is not a global minimum, but can be a local mini- 
mum. 

[18] G. Birkl, S. Kassner, and H. Walther, Nature (London) 

357, 310 (1992); I. Waki, S. Kassner, G. Birkl, and H. 

Walther, Phys. Rev. Lett. 68, 2007 (1992). 
[19] M. G. Raizen, J. M. Gilligan, J. C. Bergquist, W. M. 

Itano, and D. J. Wineland, Phys. Rev. A 45, 6493 (1992). 
[20] D. G. Enzer, M. M. Schauer, J. J. Gomez, M. S. Gulley, 

M. H. Holzscheiter, P. G. Kwiat, S. K. Lamoreaux, C. G. 

Peterson, V. D. Sandberg, D. Tupa, A. G. White, R. J. 

Hughes, and D. F. V. James, Phys. Rev. Lett. 85, 2466 

(2000). 

[21] H. Landa, S. Marcovitch, A. Retzker, M. B. Plenio, and 

B. Reznik, Phys. Rev. Lett. 104, 043004 (2010). 
[22] E. Shimshoni, G. Morigi, and S. Fishman, Phys. Rev. 

Lett 106, 010401 (2011). 
[23] G. Morigi and Sh. Fishman, Phys. Rev. Lett. 93, 170602 

(2004); Phys. Rev E 70, 066141 (2004). 
[24] Sh. Fishman, G. De Chiara, T. Calarco, and G. Morigi, 

Phys. Rev. B 77, 064111 (2008). 
[25] N. W. Ashcroft and N. D. Mermin, Solid State Physics 

(Saunders College, Philadelphia, 1976). 



